Pulsars and Gravitational Wave Detection 



George Hobbs A 

A Australia Telescope National Facility, CSIRO, PO Box 76, Epping, NSW 1710, Australia 
A Email: george.hobbs@csiro.au 



Abstract: 

The number of known millisecond pulsars has dramatically increased in the last few years. Regular 
observations of these pulsars may allow gravitational waves with frequencies ~ 10 -9 Hz to be detected. 
A "pulsar timing array" is therefore complimentary to other searches for gravitational waves using 
ground-based or space-based interferometers that are sensitive to much higher frequencies. In this 
review we describe 1) the basic methods for using an array of pulsars as a gravitational wave detector, 
2) the sources of the potentially detectable waves, 3) current limits on individual sources and a stochastic 
background and 4) the new project recently started using the Parkes radio telescope. 
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1 Introduction 

E: 



The assumption that a pulsar is a regular rotator that 
follows a predictable slow-down model forms the ba- 
sis of a powerful technique for finding its rotational 
and positional properties. This "pulsar timing" tech- 
nique (see e.g. Manchester & Taylor 1977 for a gen- 
eral review or Blandford, Narayan & Romani 1984 and 
Backer & Hellings 1986 for more details) allows the ar- 
rival times of pulses from a particular pulsar to be pre- 
dicted with great accuracy. For some pulsars that have 
spin periods of a few milliseconds (hereafter referred 
to as the "millisecond pulsars" or MSPs) the pulse ar- 
rival times can be modelled to less than a microsecond 
over many years of observations (see Figure 0. Due to 
this phenomenal precision, gravitational waves (GWs) 
with periods between days and decades should be de- 
tectable by analysing slight discrepancies between pre- 
dicted and actual pulse arrival times. As emphasised 
by Foster (1990), the effect of a GW passing a free test 
mass, such as the Earth or a pulsar, is not to move the 
mass from its coordinate position, but instead to de- 
form the space-time metric around the mass. 

Pulsar studies have already placed stringent limits 
on a GW background and have been used to rule out, 
or place limits on, some cosmic string models. As will 
be shown, pulsar studies have also placed constraints 
on postulated supermassive black hole binary systems 
in our and nearby galaxies. In section 2, we describe 
the basic framework for studying the observable effect 
of GWs on pulsar timing residuals. Section 3 contains 
a discussion on the creation and detection of a stochas- 
tic background of GWs. In section 4, we describe how 
individual sources of GWs may be detected. In section 
5 we highlight some practical issues neccessary for the 
detection of GWs. We conclude with a description 
of the current status of the Australian timing array 
project. 
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Figure 1: Timing residuals (the difference be- 
tween actual pulse arrival times and times pre- 
dicted using a spin-down model for the pulsar) for 
PSR J0437-4715. The timing residual root-mean- 
square (rms) is 400 ns. 



2 Basics 

Detweiler (1979) provided the basic framework for de- 
scribing the effect of a GW passing through the solar 
system on a pulsar's timing residuals. In brief, the 
pulsar and Earth should be considered as the ends of 
a free-mass GW antenna. In order to detect GWs the 
relative motion of the pulsar and Earth must be moni- 
tored by observing fluctuations in the pulsar's observed 
spin rate. The measured frequency v(t) of a pulsar of 
constant frequency vo, with direction cosines a, f3 and 
7 varies slightly as a GW passes the solar system as 
(Detweiler 1979) 



(a 2 -/3 2 )(hf-h P _) + 2aP(h x -h x ) 
2(1 + 7) 
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where h+' P and h^' P are the wave amplitudes in the 
two polarizations at the Earth (E) and pulsar (P). 

The observable effect of a GW is therefore to cre- 
ate pulse period fluctuations with an amplitude pro- 
portional to the gravitational wave strain evaluated at 
the Earth, h\\ x , and at the pulsar h p x . However, the 
values of h p x for widely-spaced pulsars will be uncor- 
rected, whereas the component at the Earth will be 
correlated. It is therefore possible to obtain this cor- 
related signal by combining measurements from mul- 
tiple pulsars; i.e. by using observations from a pulsar 
timing array. Hellings & Downs (1983) developed a 
simple method for determining the GW signal com- 
mon to all pulsars by cross-correlating the time deriva- 
tive of the timing residuals for multiple pulsars. We 
note that the Doppler shifts in the apparent rotational 
rates of pulsars are correlated around the sky with a 
quadrupole and higher order angular signature. More 
details about combining multiple data sets to search 
for this signature have been described by Foster (1990). 

3 Stochastic Backgrounds 
3.1 Creation 

A stochastic background of GWs can be cosmological 
(e.g. due to inflation, cosmic strings or phase transi- 
tions), or astrophysical (e.g. due to coalescing massive 
black hole binary systems that result from the mergers 
of their host galaxies). 

Models for cosmological stochastic backgrounds of 
GWs have been reviewed in Maggiore (2000). Such 
GWs can have frequencies between / ~ 10~ 18 Hz (cor- 
responding to a wavelength as large as the present 
Hubble radius of the Universe) to / ~ 10 +12 Hz (which 
corresponds to the frequency of a graviton produced 
during the Planck era and redshifted to the present 
time using the standard cosmological model). One 
mechanism for producing copious amounts of GWs 
is based upon topological defects that formed dur- 
ing phase transitions in the early universe ("cosmic 
strings"). The predicted GW spectrum due to these 
cosmic strings has an almost flat region that extends 
from / ~ 10 -8 Hz to 10 10 Hz and a peak in the re- 
gion of / ~ 10 _12 Hz. Details of the implications of 
the current pulsar timing limit are provided by Cald- 
well, Battye & Shellard (1996) who use the limit to 
exclude a range of values for the cosmic string linear 
mass density for certain values of cosmic string and 
cosmological parameters. 

Stringent limits have already been placed on the 
energy density of a stochastic background using the 
timing residuals of individual pulsars. For a flat GW 
energy spectrum that is centred on some frequency / 
and has a bandwidth also equal to / then an upper 
limit on the energy density (p g ) of a GW background 
can be obtained from (Detweiler 1979) 



Pg < 



243 n 3 f 4 



{R 2 (t)) 



(2) 



208 G 

where G is Newton's gravitational constant and (R 2 (t)) 
the rms timing residual. Romani & Taylor (1983) ob- 
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Figure 2: Predicted gravitational wave back- 
grounds. The top of the figure contains current 
and predicted limits placed on the gravitational 
wave background from millisecond pulsar timing, 
LISA and advanced LIGO. Figure provided by M. 
Kramer. 



tained an upper limits on the equivalent mass density 
in GWs in frequency ranges between 0.7 x 10 -8 and 
6 x 10 -8 Hz and concluded that, at the present epoch, 
the mass density of the universe is not dominated by 
GWs of frequency ~ 10 -8 Hz. 

The expected power spectrum for cosmological mod- 
els was shown by Blandford et al. (1984) to be 



P B (f) = 



3tt 3 / 4 



(3) 



in which p g (f) is the energy density of the stochastic 
background at frequency /, Ho = 100/irj kms -1 Mpc~ 
is the Hubble constant and il g the fractional energy 
density in GWs per logarithmic frequency interval. 
Therefore, if Q g is constant then P(f) oc /~ 5 . The 
index of this power-law (a — 5) should be contrasted 
with a — for white noise and a = 2 — » 4 expected for 
clock instabilities, ephemeris errors, interstellar prop- 
agation effects and pulsar rotational instabilities (see 
Stinebring et al. 1990). Stinebring et al. (1990) used 
seven years of observations to place rigorous upper 
bounds on the stochastic background. Using a similar 
method and seven years of data for PSR B 1855+09 
allowed Kaspi, Taylor & Ryba (1994) to place a limit 
of n gw h^ < 6 x 10" 8 . McHugh et al. (1996) pro- 
vided a more statistically sound method to obtain that 
£l gw ho < 10~ 6 . Their result is independent of the as- 
sumption of a flat spectrum for Q gm (f). 

An astrophysical background would be formed by 
GW radiation from supermassive black holes. Current 
theories suggest that galaxies contain a central black 
hole of mass 10 6 Mq . As many galaxies are observed 
to be merging, the existence of a binary black hole 
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system in a merger remnant is likely. If the binary 
loses enough energy and angular momentum then it 
may enter a regime where gravitational radiation alone 
can bring about inspiral and coalescence (Rajagopal & 
Romani 1995). During an inspiral the GWs will sweep 
through a range of frequencies. The detection of such 
events clearly depends upon the rate of occurrence and 
the amplitude and frequencies of the GWs produced. 
The possibility of detecting a stochastic background 
of such events with pulsar timing was described by 
Rajagopal & Romani (1995). 

This work was continued by Jaffe & Backer (2003) 
who find that the spectrum of a stochastic background 
of black hole binary systems has a characteristic strain 
of h c (f) ~ 10" 16 (//2/r -1 )~ 2/3 which is just below 
the detection limit from recent analyses of pulsar tim- 
ing measurements (see Figure El • Even though the 
amplitude of this spectrum was considered too high 
by Wyithe & Loeb (2003), this background is likely 
to dominate over cosmological stochastic backgrounds 
and is therefore the most likely background to be de- 
tected using a pulsar timing array. The actual back- 
ground reached can be determined from the slope of 
the spectrum which will indicate whether we are ob- 
serving a coalescing population or other sources of 
stochastic GWs. 
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Figure 3: Characteristic strain spectrum for dif- 
ferent models of a GW background formed due to 
merging black-hole binary systems. The dashed 
line gives the current best limits on the background 
from pulsar timing observations and the dotted 
line provides the expected limits from the pulsar 
timing array after ~ 8 years of operation. Figure 
obtained from Jaffe & Backer 2003. 



This characteristic strain h c (f) is defined as h c (f) = 
h(f) where is the spectral density with units of 
inverse frequency. This can be related to f2 3 (/) by 

<W) = ^CfShif) (4) 



4 Individual sources 

A periodic source of gravitational radiation will pro- 
duce a periodic shift in the pulse arrival time. Suf- 
ficient GW amplitudes and frequencies that can po- 
tentially be detected using a pulsar timing array are 
predicted to occur from a supermassive black hole bi- 
nary system in a nearby galaxy (or in the centre of our 
Galaxy). The expected signature for a supermassive 
black hole binary system with a circular orbit is given 
by (Jaffe & Backer 2003) 



h » 4.4 x 10- 17 m^ /3 P~ 2/3 d^ c — ^— (5) 



where h is an order-of-magnitude estimate of the strain 
amplitudes, m% is the total mass of the system in units 
of 1O 8 M0, P yr the observed GW period in years, d the 
distance to the emitter from the Earth in GPc and q is 
the mass ratio (q < 1). They also find that the power 
radiated along the axis of the orbit is eight times that 
for an edge-on view. The actual effect on the pulsar 
arrival times will depend upon the angle between the 
pulsar GW source and the pulsar; a pulsar lying along 
the line of sight to the GW source will experience no 
effect. 

Lommen & Backer (2001) searched for gravitational 
radiation from Sagittarius A* which had been postu- 
lated to be a massive black-hole binary system (see e.g. 
Zhao, Bower & Goss 2001). They calculated that the 
expected effect would be about ~ 10ns in the timing 
residuals of PSRs B1937+21 and J1713+0747 which 
is too small to be detectable with current data. Lom- 
men & Backer (2001) tabulated the expected timing 
residuals for postulated binary massive black holes in 
nearby galaxies assuming an equal-mass binary system 
with an orbital period of 2000 days. If we can identify 
structures with amplitudes of ~ 100 ns in the timing 
residuals then meaningful constraints can be placed 
on ~10 nearby sources. 

The expected signature for the timing residuals for 
the more general case of a coalescing, binary system in 
an eccentric orbit was presented in Jenet et al. (2003) . 
The effect depends upon the orbital parameters (in- 
cluding the orbital inclination angle), masses, source 
distance and the opening angle between the source 
and the pulsar relative to Earth. Jenet et al. (2004) 
attempted to detect GWs emitted by the proposed su- 
permassive binary black hole system in 3C66B (Sudou 
et al. 2003). The expected signature in the timing 
residuals of PSR B1855+09 are two sinusoids, one with 
an amplitude ~ 5/xs and a period of 0.88 years and the 
other of amplitude ~ 10jus with a 6.2 year period (see 
Figure 01 . The two sinusoids occur if the Earth-pulsar 
line-of-sight is perpendicular to the GW propagation 
vector as the observed timing residuals will contain in- 
formation about the source at the current epoch and 
4000 years ago (the distance to PSR B 1855+09 from 
Earth is 40001t-yr). No such signature was found. 
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Figure 4: Theoretical timing residuals induced by 
GWs from 3C66B. No such signature was observed 
in real pulsar timing residuals (Figure obtained 
from Jenet et al. 2004) 



5 Practical issues 

Pulsar timing is affected by the stability of terrestrial 
clocks, ephemeris errors and the pulsar itself. In order 
to make a definitive detection of gravity waves then a 
timing array requires a minimum of five pulsars spread 
over the sky (see Foster & Backer 1990). One pulsar is 
required to confirm the stability of terrestrial clocks, 
three widely spaced pulsars are necessary to solve for 
an error in the solar system ephemeris and the fifth 
pulsar to place an upper limit on the G W background. 
To solve completely for all the available information 
about the background, more pulsars are required. 

The pulse arrival times must be determined to high 
precision (current timing array projects are aiming to 
achieve a precision between 10 -9 and 10 -6 seconds). 
The necessary precision can be estimated from (Ra- 
jagopal & Romani 1995) 



-^-{NobsTobs 



-1/2 



(6) 



in which h is the typical strain amplitude of a de- 
tectable GW, N bs the number of arrival times mea- 
sured each year, T b a the total length of time that the 
pulsar has been observed, P w the period of the GW 
and St the typical arrival time precision. Such timing is 
achievable. For example, van Straten et al. (2001) ob- 
tained a residual root-mean-square of only 130 ns over 
40 months of observing PSR J0437-4715. With bet- 
ter instrumentation this precision should be improved. 
Measuring pulse arrival times with high precision for 
most MSPs benefits from observing at frequencies of 
~ 3 GHz or higher to counteract interstellar propa- 
gation effects (Rickett 1977). It is also essential that 
the pulsar's dispersion measure is known accurately for 
every observation. To do this, simultaneous multiple 
frequency observations are required at widely spaced 
observing frequencies. 

The pulsars chosen as part of a timing array must 
be intrinsically stable. Some pulsars show "timing 
noise" , a continuous, noise like fluctuation in the rota- 
tion rate (for example, Lyne 1999) or glitches, sudden 
increases in rotation rate (Lyne, Shemar & Graham- 
Smith 2000). However, the MSPs have been shown to 
be extremely stable over many years of observing. In 
fact, the stability of many MSPs rival that of terres- 
trial atomic time standards (for example, see Lommcn 



2002). 

The power spectrum of the time derivative of the 
pulsar timing residuals must be obtained in order to 
obtain detailed information about the stochastic back- 
ground. Methods for obtaining this power spectrum 
from the irregularly sampled pulsar timing residuals 
have been fiercely debated in the literature. Stinebring 
et al. (1990) developed a method using orthonormal 
polynomials as basis functions. This work was gen- 
eralised for multiple pulsar data sets by Kaspi, Tay- 
lor & Ryba (1994), but was criticised by Thorsett & 
Dewey (1996) who developed a technique based on the 
Neyman-Pearson test. McHugh et al. (1996) subse- 
quently showed the Neyman-Pearson test of hypothesis 
also cannot, in the general case, provide upper limits 
on an unknown parameter and suggested the use of a 
Bayesian formalism. 

When forming the power spectrum of the pulsar 
timing residuals it is important to recall that the resid- 
uals were obtained by fitting a model that includes at 
least the phase, spin-period, its derivative and posi- 
tion to the pulse arrival times. For many MSPs, fits 
have also been made for the system's orbital param- 
eters. The fitting procedure therefore removes long- 
period variations in the arrival times and hence lim- 
its on GWs are not valid for GW periods near to the 
length of the data span (see Backer & Hellings 1986). 
The transfer function of the pulsar model as a filter was 
obtained by Blandford, Narayan & Romani (1984) and 
should therefore be taken into account when studying 
relatively short data spans. 



6 The Australian pulsar tim- 



ing array 
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Figure 5: Location of the millisecond pulsars on 
the sky that are being observed as part of the Aus- 
tralian pulsar timing array project. 



Since Februrary 2004, in a collaborative effort be- 
tween the ATNF and the Swinburne and Caltech uni- 
versities, we have been observing ~20 millisecond pul- 
sars using the new 680/3100 MHz dual-frequency re- 
ceiver and a 1400 MHz receiver at the 64-m Parkes ra- 
dio telescope (see |http://www.atnf.csiro.au/research/pulsar/array| |. 
We intend to observe each pulsar at 7-10 day inter- 
vals over a period of at least five years. For many of 
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these pulsars we already have timing residual rms val- 
ues less than with observation times of 30 minutes 
or 1 hour depending upon the brightness of the pul- 
sar. For PSR J0437-4715 we currently have 166 dual- 
frequency observations spanning 125 days which give 
us, with only rudimentary processing being applied to 
the data, a timing residual rms of ~500 ns. This short 
data span already places a limit of Q gm ho < 5 x 10~ 4 
on any possible existence of a GW background. If we 
can reduce the timing residual to 100 ns over 5 years 
then the limit (from a single pulsar) will be Q gw ho < 
5 x 10~ 10 and will provide tight constraints on grav- 
itational wave backgrounds from merging black holes 
and cosmic strings. 

We have also been developing a software package, 
SUPErTempo, for processing the arrival times from 
multiple pulsars simultaneously. The current status 
of this, and other related projects, may be found on 
our web-site. 

7 Conclusion 

It is likely that gravitational waves will be detected 
within the next few years by world-wide pulsar timing 
array projects. As a by-product of these investigations 
stringent checks will also be placed on terrestrial time 
standards and the solar system ephemeris. The regu- 
lar dual-frequency observations of multiple pulsars will 
also provide valuable information about the interstel- 
lar medium. Using a pulsar array as a gravitational 
wave detector is complimentary to other searches cur- 
rently being designed that are attempting to detect 
much shorter-period GWs. 
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